PHOTOELECTRIC PHOTOMETRY OF COMET KOHOUTEK (1973f) 

Lubos Kohoutek 

1. Introduction 

Comet Kohoutek (1973f) has been observed with the 50 cm (f/15) 
reflecting telescope of the European Southern Observatory, La 
Silla, Chile, on fourteen nights between January 16 and 30, when 
the heliocentric and geocentric distances of the comet were r=0. 66 
- 1.00 A. U. and A = 0.81 - 0.96 A. U., respectively. The 40” and 
80” diaphragms were used for the photometry of the cometary head 
(EMI 6256 A photomultiplier) in the UBV system and with six inter- 
ference filters: CN 3884 R, CO"^ 4267 R. C^ 4737 R, 5172 R, 
cont. 5300 X and Na 5893 .R. 

The atmospheric conditions were good but the accuracy of our ob- 
servations was lower than usual due to large extinction (air mass 
2.5 to 4. 3) and twilight. The mean error of one measurement of 
log F in all but Na 5893 ^ filters can be estimated at +0.02, where- 
as the accuracy through the Na filter was substantially lower. 

2. UBV Observations 

Observations in UBV system are given in Table 1 (N - number of 
measurements in the respective night). They are also presented in 
Fig.l, but reduced to the distance of A = 1 A.U. and corrected to 

4 

standard circular area of D(40”) = 2. 90 x 10 km and D(80”) = 

4 

5.80 X 10 km in diameter, respectively, centered on the cometary 
nucleus. The above areas correspond to those seen in the diaphragms 
of 40” and 80” and from the distance 1 A.U. For that correction 
the mean intensity gradients in the coma were applied as deduced 
from the measurements in both diaphragms. 
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The reduced brightness of Comet 1973f in UBV system 

4 

[ A = 1 A,U. ; standard circular area of 2. 90 x 10 km (40) 

4 

and 5. 80 x 10 km (80) in diameter] as a function of log r. 
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In order to increase the time interval of observation we used the 
photoelectric measurements which were made by Vogt (1974) at La 
Silla under similar conditions. The B and V magnitudes in the dia- 
phragms of 29'* and 45”, respectively, were also corrected to 
standard circular area of 40” (A = 1 A. U. ). 

The absolute magnitudes, the photometric exponents, n, 

derived from the well-known formula 


m = m -f 2. 5 n log r 
o 


( 1 ) 


are given in Table 2 together with the respective interval of log r. 
The absolute colour indices in the 40” and 80” diaphragms differ 
from each other only slightly and within the observing errors: 

< = -''>0,(40") = <'>-=>o, (40") = < = -''>0,(80") = 

<'»- = >o, (80") = 


The comet’s brightness has been diminishing rather quickly after the 

perihelion passage, in fact more rapidly in the inner part of the head 

(n between 4. 6 and 5. 0) than in the outer part (n between 3. 6 and 

4.2). When extrapolating the best determined straight lines of V(40”) 

and B(40”) back to the perihelion, V (40”) ~ -0. 7 and B (40”) 

max max 

- -0. 5 can be obtained. 

“2 

Assuming a density law d( P ) - P ( P is the distance from the 

nucleus) for the emitting molecules and dust particles in the spheri- 
cally symmetrical head, an intensity law I( P ) ~ P ^ can be de- 
rived, which leads to the magnitude difference A m = m(40”) - 
m(80”) = 0. 75 for the diaphragms used. Our observations show 
A m > 0. 75 and continuously growing with time. The mean values 

are: A m = 0. 904, A m = 0. 817, A m = 0. 853. The difference 
V B U 

between the observed and theoretical A m- values reflects the de- 

“1 - Tf 

viation from the P law and can be expressed by P ; obviously 


53 



Table 2 

UBV Photometric Parameters 


System 


Alog r 


V(40") 

V(8C0 
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Table 3 • 

Parameters of Interference Filters 
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5f = 2 - 


( 2 ) 


2 A m 

5 log(D 2 /D^) ' 

where D , D. are the respective diameters of the diaphragms. 

^ X 

Physically it reflects (a) the production of visible radicals from 

parent molecules in this part of the inner coma, (b) the deviation 

from the regular ejection of particles from the nucleus and/or 

(c) the dissociation of the molecules by sunlight. This last effect 

can probably be neglected due to the large trajectory of molecules 

before dissociation compared with the size of the diaphragms used 

According to (2) we have: JC , 7 ^ = 0.80, =0.91, = 0 . 87. 

V x 3 U 


3. Observations through Interference Filters 

The luminosities of the cometary head in the 40" and 80" dia- 
phragms and through the interference filters (see Table 3) were 
measured using the comparison star 77 Hya, for which Hayes (1970) 
has published the energy distribution. We adopted the V magnitude 
of 77 Hya to be 4. 289 (Blanco et al. , 1968 - mean value), the V 

magnitude of the Sun -26.74 (Allen, 1973), and the flux from a star 

- 6 

of zero visual magnitude outside the atmosphere to be 3.73 x 10 
- 2 - 1 

erg. cm . s in the V region (Allen, 1973). The observed fluxes, 

F , , are given in Table 4 for CN 3884 R, C_ 4737 K and 5172 K, 

obs 2 

and continuum at 5300 X. 


Analogically to UBV data, the differences A m = m(40") - m(80") 
for molecules CN, as well as for dust in the continuum were 
found to be larger than expected according to the simple intensity 
law, and increasing with time. We found the following mean values: 

0. 77, and JC = 0. 90. The 53 QQ 

can be explained by the diminishing ejection of dust particles from 
the nucleus. Probably both effects (a) and (b) are responsible for 


^ CN, C2 


> 1 . 
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Table 4 

Observed Fluxes of Emission Bands and Continuum 
[erg. cm^ . s ] 


Filter 

Diaphragm 


Jan, 16.032 

17.027 

19.029 

21 .030 

22.031 

23.031 

24.033 

25.035 

26.031 

27.027 

28.034 
29.033 

Jan. 30.032 


CN 3884 1 
40 " 80 " 

log 

Cg 4737 ^ 

40 " 80 " 

obs 

C2 5172 X 
40 " 80" 

Cont.5300 X 
40 " 80" 

-9.234 


-9.802 


-9.290 


-10.002 


-9.209 


-9.792 


-9.305 


-9.967 


-9.154 


-9.668 


-9,21 1 


-9.913 


-9.199 




-9.250 




- 9.413 

- 9.048 

-9.901 

-9.554 

-9.444 

-9.096 

-10. 1 15 

- 9.822 


-9.067 


-9.600 


-9.122 


- 9.846 

-9.521 

-9.233 

-10.054 

- 9.814 

-9.615 

-9.357 

- 10.287 

-10.030 

-9.519 

- 9.144 
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-9.691 

- 9.602 
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-9.190 
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-9.732 
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- 10.012 
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The energy distribution of the cometary continuum has been assumed 
to be of solar type. Then, the contribution of such a cometary con- 
tinuum in the five interference filters could be calculated from the 
continuum measured in 5300 The ratio of the observed total 
flux and the calculated continuum flux in the respective filter com- 
bination is also'^given in Table 4. We find that about 95 per cent of 
the radiation transmitted by the CN 3884\R filter is the emission 
contribution and that the fraction of emissions by the 47 37 % 
and 5172 filters are 58 and 78 per cent, respectively. On the 
other hand, the observed total flux in the CO^ 4267 K filter is only 
about 4 per cent larger than that of the continuum. Our slit spectra 
of January 6 to 14 show either a missing or only very faint CO^ 
emission, so that the radiation transmitted by this filter should 
originate entirely from the cometary continuum. This fact supports 
strongly our assumption about the energy distribution of the comet- 
ary continuum. 

o 

Our measurements in the Na 5893 A filter require special care. 

The sodium doublet was very strong on the slit spectra between 
January 5 and 8 (r = 0. 33 to o.42 A. U. ) and its brightness has di- 
minished till January 12 (r = 0. 54 A. U. ) (Kohoutek, Rahe, 1974). 

The emission contribution in this filter was still about 75 per cent 
of the total observed flux on January 16 (r = 0.66 A. U.), but only 
about 18 per cent on January 17 (r = o. 68 A.U.). After January 19 
the emissions contributed only with 14 per cent assuming that the 
cometary continuum was of the solar type. If the natrium doublet 
was totally absent in r < o.7 A.U., the measurements in the 
Na 5893 .R filter would represent a colour excess of +0?^16 in the 
continuum, possibly caused by the scattering in the cometary head. 
Unfortunately, due to the very low quantum efficiency of the photo- 
multiplier in that wavelength, the measuring errors were large 
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and the above result seems to be very uncertain. For that reason 
our assumption that the cometary continuum is of the pure reflec- 
tion type has not been changed. 


The observations through the interference filters were treated in 
the same way as in case of UBV data. The observed fluxes were 
reduced to A = 1 A.U. and corrected to standard linear area in 
the cometary head, D(40”) and D(80"). First of all, the important 
ratio of the absolute brightness of the dust and of the gas coma, k = 
F , / F , could be estimated. Assuming, that the effective 
band width in the U, B, V is 680 A , 980 A and 890 A , respect- 
ively (Allen, 1973), we calculated the contribution of the cometary 
continuum to the total brightness in the respective colours. Putting 

F ~ F j and F (i) = F + F , , where F (i) is the 

o, cont o, d o o, g o, d o 

corrected flux energy in the i = V, B, U colour transformed from 
the absolute magnitude (see Table 2), we received the following 
results: 


k^(40") = 0. 69 
k^(40") - 1. 06 

D 

k^(40") = 0. 61 


k.^(80”) = 0. 66 

k^(80") = 0.84 
t> 

k^(80”) - 0. 56 


The dust coma is more concentrated toward the nucleus than the 
cometary gas. Besides, the increase of k with time in the interval 
between r 0. 6 and 1.0 A.U. was found in all colours. 


In order to know the fluxes of the emission bands, the fraction of 
their total radiation, that was passing through the respective filter, 
has to be determined. Besides, the contribution of other emissions 
in the filter was estimated, giving 7. 3 per cent in the CN 3884 R 


58 


filter, 26 per cent in the C 2 4737 R filter and 11 per cent in the 
C 5172 R filter from the total emission flux. The corrected 
fluxes of the emission bands CN(O-O) av = 0 , C^Cl-O) Av = +1 , 
and C^{0-0) Av = 0 , respectively, are plotted versus log r on 
Figs. 2, 3 and 4 and compared with the total flux measured in the 
continuum 5300 ^ (Fig. 5). The equivalent width of the 5300 R 
filter was 46 R . 

In the first approximation there is a linear dependence of log F 
on log r in both 40" and 80" circular areas for the emission 
bands as well as for the continuum. The slopes, n, and the absolute 
brightnesses ( A = r = 1 A.U.), log F^ , in the relation 

log F = log F^ - n log r - (la) 

are given in Table 5. The interval of log r was -0. 18 to -0.01 
in the 40" area and -0. 13 to 0. 0 in the 80" area. 

For continuum at 5300 R the following parameters were found: 
n(40") -4.51, log F^ (40") - -10.806; n(80") =3.81, log F^(80") = 
-10. 424. 

The photometric exponents, n , lie for the emission bands as well 
as for the continuum in the same range as in case of the UBV data 
(see Table 2). Also n(40") > n(80") can be stated with the exception 
of the C (1-0) band probably because of larger observing errors. 


4. Desorption Heats and Number of CN and C 2 Molecules 

Using Levin's theory (1943) we may express the logarithm intensity 
of the gas constituent of the coma by 
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function of heliocentric distance. See the text. 
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Table 5 

Photometric Parameters of CN and C 2 Emission Bands 


Band 

CN(O-O) 

, AV =0 

CjO-o) 

,AV =+1 

02(0-0), 

Av =0 

Diaphragm 

40 " 

80 " 

40 " 

80 " 

40 " 

80 " 

n 

4.50 

3.98 

4.62 

4.64 

5.01 

4.55 

foe F„ 2 _ 

Lerg.cm .s J 

-9.632 

-9.219 

-10.701 

-10.330 

-10.017 

- 9.614 

B 

9.51 

7.88 

9.79 

9.53 

10.67 

9.33 

f'o >2 - 1-1 

Lerg.cm .s J 

-9.650 

-9.220 

-10.721 

- 10.341 

-10.037 

-9.624 

[cal. mol” ] 

5670 

4700 

5830 

5680 

6360 

5740 

[erg . 8 J 

17.799 

18.229 

16.728 

17. 108 

17.412 

17.825 










log F = log F' - log r + B(1 - ) log e (3) 

where B = L/RT^ , L - desorption or evaporation heat, R - gas 
constant, T - temperature of the nucleus at a distance of 1 A.U. 

- Ct 

The constant a in the relation T = • r corresponds to 0.5 

if the nucleus is in a thermal equilibrium state. The constants B 

and log F^ for the CN and emission bands are also given in 

Table 5 (for ce =0.5). Both relations, (la) and (3) differ only 

very slightly from each other as shown in the example of Fig. 2. 

The desorption heat L [cal. mol was determined using T^ = 300°K. 

The mean parent molecule desorption heats of 5200 cal. mol ^ and 

5900 cal. mol ^ for the predecessor of CN and j respectively, 

differ from each other. They are larger than the value 4500 cal. mol 

which was found for CN and C„ in Comet Rudnicki (1966e) (Mayer, 

^ -1 

O'Dell, 1968), or than the value 4800 cal. mol for C„ in Comet 
Arend-Roland (1957 III) (Liller, 1965). If we use o- < 0. 5 as re- 
commended by some authors and still T = 300° K, we would re- 

° -1 

ceive B and L even larger: e.g. 6400 cal. mol (CN) and 7300 

cal. mol ^ (C_) for a =0.4, or 8500 cal. mol ^ (CN) and 9600 
-I ^ 

cal. mol (C„) for a =0.3. For that reason it seems to be more 
justified to accept = 0. 5 for Comet 1973f, assuming again T^ 
as given above. 

Following the method introduced by Wurm (1943) the number of 
molecules CN or in the cometary head can be calculated from 
the respective luminosities L: 

m 

N = L (4) 

TT fp p ( w , r) 

where f is the oscillator strength, p is the vibrational transition 
probability and p { v , r) is the solar radiation density at the given 
frequency and for the given heliocentric distance. We used the 
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molecular constants presented in Table 6 (for f see Wentink et al. , 

1964; for p see Fraser et al. , 1954) and p ( v , r) resulting from 

the solar spectral irradiance data listed by Robinson (1966). The 

absolute luminosities appearing in Table 5 were calculated from 

2 27 

the corresponding fluxes: L = 4 it F' = 2.812 x 10 F' • 

o o o 

( A = 1 A.U.). Then, the number of molecules CN and C„ which 
radiated in the cylinders extending through the comet in the line of 
sight and of diameters D(40") and D(80") were determined (Table 6). 

The ratio N (C„ 5165)/N (CN) = 1.52 or 1.46, and N (C^ 4737)/ 
o 2 o o 2 

N^(CN) = 0.99 or 0.88 for the 40" and 80" diaphragms is very 
low and comparable with the value 0.8 obtained for Comet Burnham 
(1960 II) by Arpigny (1965). Also Vanysek (1969) found for Comet 
Ikeya-Seki (1968 I) a rather low value of N(C^)/N(CN) =3.2. It is 
interesting to notice that 1960 II as well as 1973f are very pro- 
bably "new" comets and that the orbital eccentricity of 1968 I was 
also very high (0. 99915). 

There exi sts a discrepancy between N(C_) as derived from the 

\ 4737 and X 5165 bands: N (C^ 5165)/N (C., 4737) 1.6. An 

o 2 O 2 

even larger difference in N(C2) was already found for Comet Ikeya 
(1964f) by Kovar and Kovar (1965). 

5. Conclusions 

Summarizing our results based on the photoelectric observations of 
Comet 1973f we may conclude: 

After the perihelion passage the comet's luminosity was quickly 
diminishing partly due to declining ejection of both gas molecules 
and dust particles from the nucleus. The respective photometric 
parameters are given in Tables 2 and 5. Measurements in two 
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Table 6 

Number of CN and C' 


Molecules 


f D 

P (»> .r) 

N^(40") 

b 

00 

0 

0.026 0.920 

1 .90 x10”^°r"^ 

1 .74 x10^° 

4.69 x10^° 

0.0030 0.237 

5.47 x10"^°r"^ 

1 .73 x10^° 

4.14 xio^^ 

0.0030 0.731 

5.58 x10“^°r“^ 

2.65 x10^° 

6.86 xIO^^ 








diaphragms, 40" and 80" in diameter, indicate the deviation from 

the intensity law P ^ for a spherically- symmetric coma and 

- Jf 

agree with the P -law; JC decreases with increasing helio- 
centric distance, the mean 3C values lie between 0.75 and 0.91. 
The dust constituent of the coma was medium intense (k between 

0. 6 and 1. 1), and the dust particles were probably of larger size 
because of the nearly pure reflecting cometary continuum. The de- 
sorption heats for CN and C_ parent molecules were unequal 

_1 ^ 

(5200 and 5900 cal. mol ) and somewhat higher compared with 
former results. The ratio N^(C2)/N^(CN) was very low (0.9 to 

1. 5) and slightly depending on the choice of the C emission band. 
Assuming T^ = 300 K, the constant a =0.5 has been accepted 
which corresponds to the equilibrium state of the nucleus. 
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